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Abstract. The predictive power of mean- field theory is emphasized by comparing 
theory with simulations under controlled conditions. The recently developed test-field 
method is used to extract turbulent transport coefficients both in kinematic as well 
as nonlinear and quasi-kinematic cases. A striking example of the quasi-kinematic 
method is provided by magnetic buoyancy- driven flows that produce an a effect and 
turbulent diffusion. 
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1. Introduction 

What happens when fluids mix? What if a fluid is moving in a magnetized environment? 
Are there analogies between the motion of a cloud in the sky, milk in a coffee cup and the 
solar flares? The study of fluids and magnetic fields has always been a challenging branch 
of physics, leading to the development of tools of wide applicability, from meteorology 
to the study of galaxies. In particular the connection between the existence of fluids in 
motion and the amplification of magnetic fields has been investigated both analytically 
and experimentally since the beginning of the twentieth century. The generation of 
a magnetic field by dynamo action was already proposed by Larmor (1919), but a 
proper understanding of such a process requires both physical insight and a theoretical 
framework that describes the magneto-hydrodynamical (MHD) context in which the 
phenomena occur. The most common theoretical approach to MHD dynamos is the 
application of mean-field theory (Parker 1955, Steenbeck and Krause 1969, Moffatt 
1978, Parker 1979, Krause and Radler 1980). The core concept on which mean- field 
theory (hereafter MFT) rests is that turbulent systems (which include most natural 
MHD dynamos) are often amenable to a two-scale approach, where the velocity and 
magnetic fields are decomposed into mean and fluctuating components: U = U + u 
and B = B + h. The mean parts U and B generally vary slowly both in space and time, 
and capture the global, and often observable behavior of the system. The fluctuating 
fields on the other hand describe irregular, often chaotic small-scale effects. 
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Using the aforementioned decomposition the equation for the time evolution of 
the magnetic field, known as the induction equation, can be rewritten as a set of two 
equations for mean and fluctuating quantities, 

_ = Vx(c7xB)+Vx£ + 7]V^B, (1) 

— = V X [U xb) + V X {ux B) + V X (uxb)' + r]V\ (2) 

where rj is the microphysical magnetic diffusivity of the fluid (here assumed uniform), 
while £ = It X 6 is the mean electromotive force, and (u x b)' = u x b — u x b. 



Correlations such a.s u x b are at the heart of turbulent transport, and apply to 
a broad range of processes, from dynamos to the mixing of chemicals through stirring. 
Here, the key task consists in relating £ to the mean field B. Underlying symmetries 
that constrain the form of this relation are a significant help. £ is a vector so if the 
system is homogeneous and the turbulence isotropic, in what direction can it point? 
The answer is that in such a system S can have constituents pointing along the mean 
magnetic field B and the mean current density J = 'V x B/fiQ (as well as higher order 
spatial and time derivatives, see §4.2p . leading to approximations such as 

£ = aB - rjtHoJ- (3) 

The coefficients linking correlations to mean quantities are known unimaginatively as 
mean- field transport coefficients, with each one describing a distinct physical effect. In 
equation ([3]), a describes the (in)famous a effect which can drive a dynamo while rjt 
quantifies the turbulent diffusion of the mean magnetic field, and /io is the vacuum 
permeability. Note that a much more general representation of £ is given by the 
convolution integral 

£{x,t)= [ [ G{x,x',t,t')B{x',t')d^x'dt' (4) 



J to J 

with an appropriate tensorial kernel G. 

Equation ([2]) contains terms that can sometimes be neglected. Most famously, in the 

case of fluids with small magnetic Reynolds number, that is ReM = UL/rj <^ 1, or low 

Strouhal number St = Utc/L ^ 1 we can drop {uxb)' in equation ([2]) and can thus make 

an analytical calculation of the transport coefficients feasible. Under this approximation, 

known as SOCA (Second Order Correlation Approximation), equation ([2]) takes the form 

f)h 

= V X ([/ X 6) + V X (w X + r/V^b. (5) 

at 

In the limit of high Rcm (hence small St) the coefficients a and rjt reduce then to (Krause 
and Radler 1980, Radler and Rheinhardt 2007) 



T, 



a = - ju-{V xu), r]t = ju^, (6) 



where Tc is a characteristic turbulent correlation time. 
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2. The need for MFT: a status report 

2.1. Motivation 

In the astrophysical context, MFT has mainly been apphed in order to understand and 
model the origin of the Sun's magnetic field and its differential rotation (Riidiger and 
Hollerbach 2004). Direct simulations of convection in spherical shells begin to reproduce 
these phenomena to some extent (Brun et al 2004, Browning et al 2006, Brown et al 
2010, Kapyla et al 2010a), but interpreting their results properly remains difficult. This 
task is approachable only in the framework of a reasonably accurate theory. 

MFT is sometimes perceived as uncertain and even arbitrary owing to a large 
amount of parameters that are often chosen to reproduce "whatever one wants". 
Adjusting parameters at will is certainly risky and clearly not permissible in the long 
run, because it would imply a complete loss of predictive power of MFT. There are 
several reasons why the "free parameter" approach has often been adopted. Firstly, 
the conventional theory for computing turbulent transport coefficients is only accurate 
at low Reynolds numbers and not well tested at larger Reynolds numbers. Secondly, 
models of solar-like dynamos that are based on a straight-forward application of mixing- 
length ideas to computing turbulent transport coefficients (Krivodubskii 1984) do not 
reproduce the Sun: the cycle periods are too short (Kohler 1973) and the migration of 
sunspot and other magnetic activity is poleward, not equatorward, which is also found 
in direct numerical simulations (Gilman 1983, Kapyla et al 2010a). 

In this situation it is sensible to reduce ones ambitions and to focus on phenomena 
that are seen in direct simulations of simplified systems which are nevertheless relevant 
for understanding the Sun. A useful goal consists then in reproducing such phenomena 
by mean-field models, thus obtaining a chance to trace down the reasons for discrepancies 
between both representations. This will be exemplified in ^ First, however, we shall 
summarize the basic saturation phenomenology of mean-field dynamos. 

2.2. Saturated dynamos and magnetic helicity fluxes 

A recent discovery that is now well explained by MFT is the slow saturation behavior 
of an dynamo in a triply periodic box (Brandenburg 2001). Such systems are 
unphysical, but they make good test problems due to the ease of capturing them both 
numerically and analytically. For early times, both the mean and the fiuctuating fields 
grow exponentially - as expected from kinematic theory. However, when the small-scale 
field becomes comparable to the equipartition value, i.e. ~ B^^, the growth changes 
its nature: the fiuctuating field saturates while the mean field, well below equipartition, 
continues to grow albeit extremely slowly. Finally, after multiple microphysical resistive 
times, the mean field itself reaches a steady, super-equipartition state; see Fig. 9.4 of 
Brandenburg and Subramanian 2005). 

This behavior is one aspect of the "catastrophic" a-quenching phenomenon, and 
has come to be understood in terms of the magnetic helicity density, h = A ■ B, 
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and the magnetic a effect of Pouquet et al (1976), where the growing magnetic a, 
ajvi = Tj ■ b/3p, is subtracted from the a of equation (|6]), now marked (kinetic), 
so that the net a = ax + would be reduced and the dynamo growth halted. This 
result has been extended to the dynamic a quenching phenomenology (Kleeorin and 
Ruzmaikin 1982, Field and Blackman 2002, Blackman and Brandenburg 2002), where 
the mean magnetic helicity in the fluctuating fields, /if = a ■ 6, is used as a proxy for 
the current helicity, j ■ 6 ~ k'^hf, with fcf being the wavenumber of the energy-carrying 
eddies. In a system that disallows the divergence of the magnetic helicity flux, such 
as a triply periodic domain, the time evolution of hf and the resulting dynamical a 
quenching equation can be written as: 

^ = -2e.B- 2r^j . b, (7) 

da n i2 f - VtJ ■ B a- orA /„v 
* = -2"*"' [ Bl + ^ j ■ 

The three phases of the dynamo in a triply periodic domain can now be understood. 
First the fields grow exponentially and the magnetic a effect grows with them until 
au ~ — Q^K- This occurs rapidly enough that still /im ~ —hf and so B / B^q — ^/ki/h < 1 
(Brandenburg 2001), where ki is the smallest possible wavenumber in the domain. 
During the resistive phase, the fluctuating flelds are nearly steady but there is still 
a small excess of a over r]ki needed to replenish the fleld in the face of resistivity. This 
phase ends only when the time evolution of the total magnetic helicity reaches a steady 
state, which occurs when {J ■ B) = 0, or for B / B^q — \fkilk'\ > 1. 

A short exponential growth phase yielding only weak mean flelds poses severe 
problems for astrophysics as the subsequent resistive growth phase is generally 
prohibitively long. Real systems however allow for fluxes of magnetic helicity across 
their borders, and/or show spatial variations in the a effect, particularly regions where 
a has opposite signs separated, say, by an equator. This raises the possibility that 
the magnetic a will be exported from the system or transported to the equator and 
destroyed. Research into such transport is recent, but has already shown conclusively 
that there is a flux of hf and that a larger residual a effect results due to it. It is not 
yet clear how large an effect this has on the flnal mean fleld strength. 

3. Predictions versus realizations 

In this section we contrast the results of some computer realizations with corresponding 
mean-fleld predictions. We discuss examples from both linear and nonlinear regimes. 

3.1. Parity and dependence on boundary conditions 

A relatively old example is the emergence of oscillatory dynamo solutions in local 
models of accretion discs (Brandenburg et al 1995). Here, turbulence is driven by the 
magneto-rotational instability which generates a negative a effect in the upper half of the 
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disc (Brandenburg et al 1995, Ziegler and Riidiger 2000, Brandenburg 2005a, Gressel 
2010). According to MFT, this negative a, together with pseudo- vacuum boundary 
conditions (that is, normal field conditions), predicts travehng wave solutions that are 
symmetric about the midplane and migrating toward the boundaries (Brandenburg 
and Campbell 1997). Conversely, when the boundary condition is changed to a perfect 
conductor condition, one expects non-oscillatory solutions that are antisymmetric about 
the midplane. Again, this dependence is born out by simulations (Brandenburg 1998). 

3.2. Onset of convective dynamo action depending on rotation rate 

Large-scale dynamos due to turbulent convection are of particular interest in 
astrophysics. According to MFT, rotating inhomogeneous (usually due to stratification) 
turbulence leads to the generation of kinetic helicity and thus an a effect which should 
enable the generation of large-scale fields. However, numerical simulations of rotating 
convection at first failed to show large-scale dynamo action (e.g. Brandenburg et al 
1996, Cattaneo and Hughes 2006). Erroneously low values of a determined with what 
is now often called the imposed-field method, where a uniform magnetic field Bq is 
applied and one determines, from the simulation, the mean electromotive force, {ux b), 
where mean fields are defined as volume averages, (e.g. Cattaneo and Hughes 2006; see, 
however, Kapyla et al 2010b) seemed to confirm that the a effect does not work. On the 
other hand, when computing turbulent transport coefficients for convection using the 
test-field method (see §4.ip it was discovered that, as the rotation rate in non-shearing 
runs increases, the a effect increases and turbulent diffusion, rjt, decreases (Kapyla et al 
2009a). Mean-field models (hereafter MFM), using these properly determined transport 
coefficients then suggested that a large-scale dynamo should be excited when the Coriolis 
number, defined as the ratio of the rotation period and the convective turnover time 
exceeds a value of ~ 4. Subsequently, direct simulations in this parameter range just 
confirmed the existence of a large-scale dynamo (Kapyla et al 2009b). Again, this 
demonstrates that already kinematic MFT has predictive power and that very likely 
MFMs can give useful and new information even about more complex systems. 

3.3. Helicity considerations 

The helicity considerations outlined in §2.2l provide further predictive power. The shear- 
induced (non-diffusive) magnetic helicity ffux, introduced by Vishniac and Cho (2001), 
has been particularly important in explaining the existence or the absence of a large- 
scale dynamo. For example, Tobias et al (2008) presented simulations of convection 
with vertical shear where no large-scale dynamo was excited although the necessary 
ingredients (inhomogeneity, rotation and shear) were all present. However, in that 
case the shear-driven magnetic helicity ffux is directed along the periodic x direction 
and no net ffux out of the system could occur. Using instead, in an otherwise similar 
setup, normal field boundary conditions, which do allow a net fiux, Kapyla et al (2008) 
showed that a large-scale dynamo does exist and indeed saturates at near-equipartition 
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Figure 1. Left: Energy in the horizontally averaged streamwise magnetic field from 
two convection simulations with vertical shear Uy{z) and either normal field (solid line) 
or periodic (dotted line) magnetic boundary conditions in the x direction. Adapted 
from Kapyla et al (2008). Right: Evolution of the ratio of the energy contained in 
the large-scale field to the total magnetic energy for open (vertical field) and closed 
(perfect conductor) boundaries. Note that large-scale dynamo action is only possible 
with open boundaries. Adapted from Brandenburg (2005b). 



field strengths; see Fig. [H where we also show the effect of open versus closed boundaries 
for forced turbulence (Brandenburg 2005b). 

Another issue approachable through magnetic helicity considerations is the 
convergence problem of the MRI, i.e. the steep decrease of the turbulence level at low 
magnetic Prandtl numbers, Pr^, in fully periodic setups (e.g. Fromang et al 2007). 
An otherwise similar setup, however, that does allow a magnetic helicity flux through 
the vertical boundaries produces indeed strong large-scale dynamo action, roughly 
independent of the value of Prjv/ (Kapyla and Korpi 2010). 

4. Computing mean-field transport coefficients 

In view of such success stories there should be an unbroken interest in MFMs both 
because of their descriptive capabilities and their predictive potential, but we have 
to realize that there are serious shortcomings of MFT that have persuaded many 
researchers to (re)turn to global direct numerical simulations instead of designing 
improved MFMs. This critical stage of MFT can be characterized by the following 
observations: 

• The limitations of analytic approaches to calculating mean-field coefficients are 
clearly too restrictive as the interest has moved from pointing out the qualitative 
existence of certain effects to the quantitative reproduction and prediction of such 
processes. This is due to 

— the obvious insufficiency of SOCA in astrophysical contexts as usually Rcm 3> 
1 and St ^ 1 

— the unclear aspects of closure approaches like the r approximation (Radler and 
Rheinhardt 2007) 
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— the need of knowledge of velocity correlators like UiUj . . .Un even in mathemat- 
ically well established (systematic) higher-order correlation approximations 

• It is obvious that MFMs for reahstic setups with predictive abilities need to employ 
transport coefficients that are 

— fully tensorial 

— position dependent 

— dependent on the mean quantities themselves, i.e. nonlinear 

— non-local and non- instantaneous 

— including magnetic background fluctuations 

There is no longer any chance for obtaining powerful models by employing a few 
scalar coefficients, the basic structure of which can be derived analytically leaving 
a few free parameters to be adapted properly. 

To flnd a way out of this impasse one might look at how in other fields of 
physics/engineering, modeling and simulation of rather complex systems are being made 
possible if the full resolution of the microphysics is not affordable. Let us choose as an 
example the mechanics of elastic materials, say metals. Their elastic properties can in 
principle be derived from the microphysics of their lattices, but it will perhaps never 
be possible to simulate the static and dynamic behavior of, say, a bridge by solving 
quantum physical lattice equations. Instead, one relies upon the equations of continuum 
mechanics in which the lattice physics enters via macrophysical material properties 
like Young's modulus and Poisson number (sufficient for an isotropic material). These 
are typically obtained by measurements in a series of standardized experiments with 
test bodies having simple geometries and being subject to clearly defined boundary 
conditions. Of course, for such an approach to be successful a certain locality of 
the microphysics processes is necessary, which can be expressed by the principle that 
neighboring material elements of a structure "communicate" only via forces at their 
common borders. 

A widely analogous procedure with respect to the task of calculating transport 
coefficients for a certain type of turbulence would consist in creating it in a (small) 
test volume with well defined boundary/environmental conditions (like a penetrating 
magnetic and/or gravitational field) and to determine then the wanted coefficients 
somehow by measurements. Then a major theoretical challenge consists in specifying 
the set of experiments needed to find just these coefficients and in prescribing the 
computational recipe for extracting them from the measured quantities like fluctuating 
magnetic flelds and/or velocities. 

This program has been implemented by the so-called test-field methods (Schrinner 
et al 2005, 2007) with the modiflcation that the physical experiments are replaced 
by numerical ones and that occasionally, in some of them, the ruling physics has 
been modifled from the canonical one. Clearly, the most important difference to the 
continuum mechanics scheme lies in the fact that just the same equations whose direct 
simulation was felt to be not affordable, and that just created the need for a MFM, have 
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now nevertheless to be simulated within the numerical experiments. However, in two 
aspects the test-field approach can still be advantageous. Firstly, the "experimental" 
volumes can represent small sections of the object which is to be globally analyzed. 
Hence, much finer structures can be resolved for the same numerical effort. Secondly, if 
a MFM is established once, it can thereafter be utilized for long-term global simulations 
which would otherwise be prohibitively expensive. 

The aforementioned locality of the actual physics has here to be required with 
respect to correlation properties of the underlying fluctuating flelds, say a turbulent 
flow. For all conceivable astrophysical situations this condition can hardly be considered 
too restrictive. In practice, correlation lengths and times are the relevant quantities to 
be considered in deflning the simulation box size and the integration time. 

Computing turbulent transport tensors like a and r/ (see Eq. ([9]) below) has now 
been done with variable success over the last 20 years. Utilizing the imposed-fleld 
method for ck has either led to the conflrmation of well-known results (for example a 
positive horizontal a effect in the upper layers of convection in the Northern hemisphere), 
or to the prediction of as yet unknown results (e.g. a reversed sign of the vertical a 
under the same conditions; see Brandenburg et al 1990) later conflrmed by theoretical 
calculations (Ferriere 1992, Riidiger and Kitchatinov 1993). 

After having explained the new test-fleld method in the next section, particular 
applications considering the non-locality of turbulent transport in space and time will 
be discussed in §4.21 



Test-field method 

Let us return to equation (jj]) for the fluctuating magnetic fleld. The wanted mean 
electromotive force £ = u x b is obviously a linear and homogeneous functional of B 
and we may therefore write the ansatz 

S = aB- rjVB, (9) 

strictly valid for stationary mean flelds depending only linearly on position. The 
components of o: and r] can be found by the following procedure: 

(i) solve 




- 7]V'^b^ - curl \Uxb'' + {ux b'')'] = curl (n x B'') 



for given u, U and test fields B'^, k = 1, . . . , N, 

(ii) calculate S'^ = uxb^, 

(iii) determine the components of ck, r] by inverting 

S'' = (xB^ - rjVW. (10) 

The solution is unique, if chosen appropriately and the test flelds B^ are sufficiently 
independent. Since we "look at" the given velocity fields u, U not only from one 
perspective like in the case of the imposed-field method, but obtain instead different 
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views represented by the different test solutions b^, the test- field approach could 
be characterized as "holographic" instead of "photographic". Indeed, the whole 
information needed to specify an ansatz like ffTOl) is extracted from the velocity fields. 

For a number of explicitly given flows like those introduced by Roberts (1970) and 
Galloway and Proctor (1992), exact agreement of the determined tensors with SOCA 
results. In the case of the Roberts flow there is agreement even with an analytic result for 
arbitrary magnetic Reynolds numbers (Radler et al 2002, Rheinhardt and Brandenburg 
2010, Radler and Brandenburg 2009). 

The method has now been applied to a number of different flows ranging from 
homogeneous forced turbulence without shear (Sur et al 2008, Brandenburg et al 2008a) 
to cases with shear (Brandenburg et al 2008b, Mitra et al 2009) and to inhomogeneous 
turbulence in stratified disks (Brandenburg 2005a, Gressel et al 2008) as well as 
convection (Kapyla et al 2009a). 

4-2. Nonlocality in space and time 

Because the test fields are not native to the system, they can disentangle effects that 
can not otherwise be distinguished. However, for the same reason they can introduce 
temporal or spatial scales that are again not native to the system. 

With respect to temporal scales the consequences of this mismatch can be seen in 
the memory effect: Consider a dynamo system with a growing mean field. Turbulence 
creates a fluctuating field from the mean one, and as long as it survives it contributes to 
the electromotive force. If the mean field is growing, the small scale field created in the 
past is weaker than it would be when created in the present. Thus, if the time behavior 
of the test fields is different from that of the "real" mean field, the ex. and rj tensors from 
the test-field method will not be the actual ones which rule the evolution of the mean 
field. A similar problem occurs when the spatial scales of the test-fields do not coincide 
with the spatial scale of the mean field to be modeled, due to non-locality in space. 

When the proper scales of the system are known, corresponding test-fields can 
be used. Otherwise, the scales of the test fields, say, wavevector k and frequency u 
are considered independent parameters and the test-field method provides OL{k,u) and 
r]{k,uj) which exhaustively describe the response kernel introduced in (jlj). 

The memory effect is demonstrated in Fig. |2] for the case of the Roberts flow (for 
details, see Hubbard and Brandenburg 2009). In the left plot, we see the difference 
between the growth rate of a dynamo to that calculated from the dispersion relation 
using cx and t} determined by the test-field method with steady test fields. We can 
reconcile these growth rates by deriving them all from a proper kernel which can be 
established by employing a set of test fields with different time dependencies. 

The memory effect and non-locality in space have been studied using the integral 
kernel technique in Hubbard and Brandenburg (2009) and Brandenburg et al (2008a), 
respectively. Using test-fields that oscillate sinusoidally in time, the Fourier transforms 
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Figure 2. Left: Rbm dependence of the dynamo growth rate for the Roberts flow as 
obtained from a direct calculation (Agrowth) compared with the result of the dispersion 
relation, Adisp = \ctk\ — {rj + rit)k'^ , using a cube of size L^; ki = 2Tr/L, kf = \f2k\. 
For this range of Rcm, the most unstable mode has the largest possible wavelength 
(A; = fci). Right: Real and imaginary parts of ol(lS) for fc = using the Otani (1993) 
MW+ flow with Rejv/ = 1. Normalization given by — Uq. Inset: scaling of Ima 
near the origin with slope 0.2, in agreement with the results of Hughes and Proctor 
(2010). Adapted from Hubbard and Brandenburg (2009). 



a{u)) of the turbulent transport integral kernels a{t) were found to fit the form: 

ajuj) _ 1 - iujTg ^^^^ 
ao {1 - iuT^y + u^T^' 

where Tq, is the memory time of the flow and Ua is a fit parameter of order r~^. In 
turbulence is comparable to the turn-over time but in steady flows it can approach 
microphysical resistive time scales. In the right panel of Fig. |2]we present such a fit for 
the MW+ flow of Otani (1993), with an extra spike at w/wf = 2, a natural frequency of 
this flow. The slope of the imaginary part at the origin, dim d;/dco'|^=0) represents the 
coefficient of the first order term with respect to an expansion in time. Its value of 0.2 
is in agreement with that found by Hughes and Proctor (2010). The real-space integral 
kernel corresponding to (fTTj) reads 

a(r) = ao6(r)e"^/^" cosw^r, (12) 

where r is the time distance to the instant of consideration, and B is the Heaviside step 
function that preserves causality as the time integral kernel must only consider the past. 

The Fourier transform of the spatial integral kernels is somewhat simpler, being fit 
by a Lorentzian: 

^= and (13) 

ao l + [aak/kiY r]to 1 + (a^fc/fcf)^ 

whose amplitude is nearly independent of Rcm for Rcm ^ 1 and the k dependence is 
reasonably close to quadratic for k/kf < 2; see Fig. |3l Here, Oq, ~ 2a^ ~ 1 are coefficients 
of order unity. The corresponding spatial integral kernels are simple exponential decays: 

^^^^^^^-ikMC ^^(^) = !^e-feKK, (14) 




Figure 3. Left: Dependence of the normalized values of a (dashed or red line) and 
77t (soHd line) on ReM for fc/fcf = 0.2 and Re = 2.2. Adapted from Sur et al (2008). 
Right: Dependences of the normalized a (dashed or red line, small symbols) and fjt 
(solid line, bigger symbols) on the normalized wavenmxiber fc/fcf for turbulence forced 
with fcf/fci = 5, Rcm = 10 (squares) and fcf/fci = 10, Rej\/ — 3.5 (triangles). Lines 
give the Lorentzian fits (|T3| . Adapted from Brandenburg et al (2008a). 



where C, is the measure of distance from the point of consideration. 
5. Prom linear to nonlinear 

As outlined in §4.H a mathematically well-founded procedure is available by which the 
tensors ex. and r] can be obtained when the velocity is given. Naturally, the question 
arises of how one should proceed when the mean magnetic field has already acted upon 
this velocity. Inspecting equation (|2]) it can be concluded that considered as a 
functional of U and -B, is always linear and homogeneous in the latter, irrespective of 
the effects U was subjected to, and whether a mean field had already acted upon it 
or not. That is, the presented method continues to be valid without modification and 
as a tribute to this extension it is called the quasi-kinematic method. The turbulent 
transport coefficients are of course now depending on B, but this dependence is entirely 
conveyed by the dependence of on 

Clearly, a dynamically effective mean field now represents an additional preferred 
direction. As a consequence, for an isotropic hydrodynamic background and a uniform 
B, the formerly isotropic tensor at adopts now the shape 

aij = aiSij + a2BiBj, i,j = 1, 2, (15) 

with B the unit vector in the direction of B. If this is, say, the x direction we 

get ail = ai + 0(2 and 0:22 = ai- Both coefficients are of course functions of B. 

Since £ = auB, the effective scalar a effect is just given by an. As an example, 

— 4 

the a quenching characteristic for the Roberts flow was determined exhibiting a B 
asymptotic dependency, cf. Rheinhardt and Brandenburg (2010) and Fig. HI This result 
is at odds with theoretical predictions, although it agrees with data for a forced ABC 
flow (Sur et al 2007). 
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Figure 4. Left: Variation of a witli B for the driven Roberts flow with Rej\/ = 
l/2\/2 w 0.35 and Ptm — 1- Adapted from Rheinhardt and Brandenburg (2010). 
Right: ReM-dependence of the normalized a and fjt in the saturated state with 
B « Bcq- Adapted from Brandenburg et al (2008c). 

Things become more involved if the direction of the mean current density J 
enters as a second preferred direction. A situation in which this comphcation is 
circumvented, without being as simple as the former one, is given by the dynamo 
due to homogeneous isotropic helical (forced) turbulence. Here, the dynamo solution 
is a Beltrami field with B \\ J and constant modulus. Hence, J is not providing 
an additional preferred direction and the form of equation ( fTSi) remains valid. The 
solution has always Beltrami shape, regardless at what level it eventually saturates. 
Consequently, a and ?7t are independent of position for any field strength. The growth 
rate of the dynamo is given by A = \ak\ — [rj + rit)k^ and should approach zero in 
the course of saturation. Under these conditions it is possible to confirm the quasi- 
kinematic method in the way that a{B) and rjtiB) are determined in the saturated 
stage and checked for consistency against A = 0. Indeed, this could be demonstrated 
to high accuracy for different values of Re^; see Brandenburg et al (2008c). Figure H] 
shows a{B) and rit{B) as functions of Rcm in the saturated state with B fa B^q. 

6. Quasi-kinematic method for magnetic buoyancy-driven flows 

So far we have been dealing with situations in which a hydrodynamic background was 
provided independently and the mean magnetic field occurred as an additional, at most 
coequal participant. But what about cases in which the turbulence itself is a consequence 
of B7 Examples are the magneto-rotational instability and the magnetic buoyancy 
instability (see below). Clearly, those setups do not know a kinematic stage on which 
the influence of B is negligible. One might worry that in such a situation the quasi- 
kinematic test-field procedure fails (Courvoisier et al 2010). However, equation ([2]) 
continues to be valid and hence all conclusions drawn from it. Consequently, the quasi- 
kinematic method should be applicable. The only peculiarity occurring here is the fact 
that all components of a and rj vanish for < i? < ^threshold; because a fluctuating 
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velocity (and magnetic field) develops only after the instabilities have set in. 

Let us now consider the magnetic buoyancy instability. It has been proposed by 
Moffatt (1978) that, once the dynamo-generated magnetic field in the overshoot layer 
of the Sun reaches appreciable strengths, this instability can set in and govern the 
dynamics thereafter. The buoyancy instability of a localized flux layer in the presence 
of stratification and rotation was later studied in detail by Schmitt (1984, 1985). A 
necessary but not sufficient condition for this instability is 



which essentially means that the magnetic field decreases faster with height than density. 
Using the imposed-field method, Brandenburg and Schmitt (1998) performed numerical 
calculations to determine the a effect of the resulting turbulence. Here we determine 
all components of a and 77 using the quasi-kinematic test-field method wherein mean 
fields are defined as xy averages. 

Our setup is similar to that described in Brandenburg and Schmitt (1998). The 
computational domain is a cuboid of size — 1 < x < 1, — 3 < ?/ < 3, —0.5 < z < 1.5, 
with gravity pointing in the negative z direction, and rotation ft making an angle 6 
with the vertical. The pressure scale height is Hp = 1 (half height of the box). The 
base state is a polytrope with index m = 3 (the adiabatic value is here 3/2), so that it 
is stable to convection. The initial condition comprises a horizontal magnetic layer of 
thickness = 0.1 with the profile. 



where the ratio of the Alfven speed to the sound speed is wao/cso = 0.5. We modify 
the base state such that the density profile remains polytropic but the entropy profile 
is adjusted to obey magnetostatic equilibrium. The initial velocity consists of about 20 
localized eddies with Mach numbers of about 10~^ at 2 = 0.1 in the xy plane. We use 
stress-free boundary conditions for the velocity and the vertical field condition for the 
magnetic field, whereas with respect to entropy we keep the temperature at the top and 
the (radiative) heat flux at the bottom constant. All calculations have been done with 
Pr = PrAf = 4 on a 64^ grid. Fig. [Sti gives the time series of the volume averages (m^) 
and (B^) and Fig. [5]d the evolution of the mean field B^. There is a short exponential 
growth phase followed by decay on a resistive time scale. 

When it comes to applying the test-field method, an aspect not discussed up to now 
is the intrinsic inhomogeneity of the flow both due to stratification and the background 
magnetic field itself. Within kinematics, that is without the background field, no specific 
complication is connected to this as from the stationary version of equation (jl]) a and rj 
emerge straightforwardly in a shape expressing inhomogeneity, that is, a{x, x'), r]{x, x') 
or, equivalently, a{x,x — x'), r]{x,x — x'). When subjecting the latter to a Fourier 
transform with respect to their second argument, we arrive at a.{x,k) and f]{x,k). In 
our case, harmonic test fields with different wavenumbers k in the z direction can be 
employed to obtain cx{z,k) and •f}{z, k). 




(16) 




(17) 
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Figure 5. (a) Evolution of volume averaged velocity and Bx squares, (b) A 'butterfly 
diagram' for B^- Panels (c) through (f) show selected tensorial components of a and t] 
for different wavenumbers k as explained in the legend in (d). Plots have been obtained 
by time averaging over the interval indicated by the dashed-dotted lines in (a). 



In the nonlinear situation, the Green's function approach remains vahd if £ is 
considered as a functional of U and B which is then linear and homogeneous in the 
latter (cf. However, we have to label G by the B actually acting upon U, that is, 
G{x,x'; B), and can thus only make statements about the transport tensors for just 
the particular B at hand. Hence, the tensors have to be labelled likewise: ct{z, k; B), 
f]{z, k; B). As our initial mean magnetic field is in the y direction, the instability will 
generate a and we are mainly interested in the coefficients Oy^, cxyy, rjy^ and rjyy with 
rank-2 tensor components rjij = —rjiksejks] they are shown in Fig. [5l (c) through (f). 

Our goal is now to confirm that the relationship between £ and B taken directly 
from the DNS can be represented by equation (ITOj) with the transport tensors deter- 
mined using the quasi-kinematic test-field method. In mathematical terms 

£i{z; B) = Rej ^ [&ijiz, k; B) - ikr]ij3{z, k; B)] cf^ exp{iTckz)^ (18) 

k 

with 2c^^^ = J B{z) exp{—iTckz)dz and = signifying the question whether both sides are 
indeed equal. We find that we can reasonably reconstruct the mean emf by truncating 
the infinite Fourier series already at A; = 8. The result of the assembly of the £ as 
formulated on the right hand side of (|T8l) is presented in Fig. |6fc and turns out to 
be a faithful reproduction of £ shown in Fig. [6^ especially during the exponential 
growth phase. We conclude that the quasi-kinematic test-field method may be used 
for correctly calculating transport coefficients even in the presence of inhomogeneous 
turbulence driven by an initial mean magnetic field. 
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Figure 6. (a) The mean emf £xit,z) calculated from the horizontal average u x b. 
(b) Reconstruction of £x{t, z) using only the fc = contributions in (fTS)) . (c) Same as 
(b) but using all contributions fc = 0,l,...,8. 



7. Conclusions 

Mean-field theory has still a lot to offer in terms of new effects and quantitative precision 
by combining analytics with numerics in parameter regimes that were previously 
inaccessible. The list of examples goes on and on; here we have only mentioned some 
of the most striking cases. The unmentioned ones concern, e.g., the Reynolds and 
Maxwell stresses that have important applications in accretion disks (Blackman 2010) 
and possibly sunspot formation (Brandenburg et al 2010). One may anticipate that all 
these aspects of mean-field theory will soon gain in significance in our voyage toward 
understanding astrophysical dynamos. 
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